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Abstract 

An approach to optimizing the energy efficiency of a Heating, Ventilating, 
and Air Conditioning (HVAC) system is presented that utilizes computa¬ 
tional predictions of the effect of heat load distribution on moist air temper¬ 
ature, density, and humidity variation. Lumped-HVAC (L-HVAC) is a new 
lumped parameter code that couples fluid transport, energy transport, ther¬ 
modynamics, and psychrometrics in an HVAC system. This code contains 
a nonlinear implicit solution algorithm for steady-state and transient calcu¬ 
lations for flow resistance, water mass balance, and energy conservation. L- 
HVAC has been validated using a simplihed analytical model, the commercial 
lumped parameter code SINDA/FLUINT, and experimental measurements. 
Steady-state calculations for a single-room system suggest an order of mag¬ 
nitude greater energy savings using a variable chiller power control approach 
compared to control damper and variable-drive fan approaches. L-HVAC was 
also applied to predict that the fraction of latent to total heat load influences 
the steady-state system temperature by up to 0.4°C for the example system 
in this study. 

Keywords: lumped analysis, energy management, HVAC system, 
psychrometrics 


1. Introduction 

Heating, Ventilating, and Air Conditioning (HVAC) systems comprise a 
large portion of energy use in today’s buildings. In 2001, the Department of 
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Energy [1] stated that 31.2% of energy use by U. S. residential buildings is by 
HVAC equipment, which equates to 355.7 billion kWh. In some specialized 
cases, HVAC energy use is more signihcant: cooling energy in data centers 
now surpasses the computing energy! [2, 3] 

The traditional HVAC design approach involves static adjustment of out¬ 
let dampers to match required airflows for assumed static loads. This ap¬ 
proach provides a reasonable approximation of required airflows for cases of 
hxed loads, but in many instances the loads fluctuate greatly with time. Fluc¬ 
tuating transient loads are apparent in classrooms, training centers, kitchens, 
and eating areas. These transient loads result in periods of overcooling due 
to small loads followed by insufficient cooling due to large loads. In addition, 
multiple solutions exist for providing the adequate amount of airflow into the 
zones, but these result in different amounts of energy consumption. 

Various approaches have been proposed or implemented in today’s build¬ 
ings to help save energy use in transient loading. The typical office building 
HVAC system applies a Variable-Air-Volume (VAV) approach to reduce cool¬ 
ing airflow into a room via thermostat control feedback [4]. This approach 
solves the problems associated with static damper placement and is inexpen¬ 
sive, so it is commonly implemented in newer office buildings. 

The insufficiency of VAV systems lies in their limited implementation to a 
single building zone. The HVAC system is actually a complex arrangement of 
fluid systems where VAV adjustment for a single zone affects the airflow into 
other zones. This interdependence between each of the control devices in a 
typical HVAC system results in a complex control scheme that is difficult to 
optimize for minimal energy usage. As a result, several studies have employed 
optimization techniques to HVAC system models such as genetic algorithms 
[5] and fuzzy interference [6, 7]. 

This study provides an alternative approach for energy minimization in an 
HVAC system that is straightforward to implement in most HVAC systems. 
The general methodology to this approach is as follows: 

• An inhnite number of HVAC designs provide adequate cooling to all 
zones in a system with hxed heat loads. The amount of energy used 
by the system depends on the chosen design. 

• Steady-state simulations of an HVAC system allow for predictions of 
HVAC system energy usage for a large number of unique designs. 

• An optimization routine pinpoints the design that provides adequate 


2 



cooling yet minimizes energy usage for a system with known fixed heat 
loads. 

• Knowledge of the instantaneous loading state in a transient system 
allows for estimation of the control conhguration that minimizes energy 
for that state using multidimensional interpolation. 

This paper focuses on the first two bullet points above in describing a simple 
predictive code that calculates the energy use and outcome of an HVAC de¬ 
sign for a user-specihed system under steady loading. Subsequent studies will 
apply an optimization routine to predict the HVAC design featuring minimal 
energy use, and then translate the optimal steady-state control conhgurations 
to transient systems. Note that this study aims to minimize energy, while 
other studies use electricity cost; this extension is straightforward provided 
electricity rates are available. Figures 1 and 2 show that the optimization 
routine hrst creates a response surface by selective sampling of control pa¬ 
rameters, and then it optimizes the the control parameters using the response 
surface and a known optimization method. This approach, described in many 
texts (e.g. [8]) and applied in numerous studies outside of the HVAC field 
(e.g. [9, 10]), is similar to that by Fong et ah [11, 12], where their study 
optimized the parameters involved in a specihc HVAC application; and Sny¬ 
der and Newell [13], who calculated the effects of building mass in a lumped 
fashion. This study provides a more general approach to modeling of HVAC 
systems, and therefore the code discussed here is intended to be used as a 
tool that can model most practical HVAC systems of interest. 

This paper focuses on the HVAC modeling code instead of the optimiza¬ 
tion algorithm. To be useful, the modeling code requires reasonable tempera¬ 
ture and humidity predictions while still being computationally inexpensive. 
The code also needs to be adaptable to most HVAC systems of interest. 
Therefore, a lumped parameter code is proposed as a good balance between 
accuracy and computational efficiency. This code, called Lumped-HVAC (L- 
HVAC), provides reasonable predictions of transient and steady-state moist 
air properties in an HVAC system with known equipment and loads as will 
be demonstrated later in this report. 

2. Code Algorithm 

L-HVAC divides the moist air in the entire HVAC system into lumps of 
moist air mass. These lumps are dehned as a fixed control volume within the 
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Figure 1: General algorithm for response surface generation. 
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Figure 2: General algorithm for system optimization. 
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system, and at any given time the moist air in the control volume has uniform 
thermodynamic and psychrometric properties. Moist air flows between these 
lumps using dehned paths that represent a tube of known cross-sectional size, 
length, and minor loss factor. 

The approach to investigating a fluid system in this manner is not unique 
to L-HVAC: the code SINDA/FLUINT [14] applies the same formulation to 
a variety of fluid systems, and it is especially useful for systems where heat 
loss from the fluid to the surrounding environment is vital. The difference 
between L-HVAC and SINDA/FLUINT lies in L-HVAC’s coupled water bal¬ 
ance, psychrometric, and energy exchange calculations. SINDA/FLUINT 
contains external subroutine input capability, so it is possible that the com¬ 
mercial code could be adapted to allow for the same calculations presented 
here. In this study, SINDA/FLUINT is used in this study for validation of 
simple systems. 

It should also be noted that the EnergyPlus [15] software package is widely 
used to model the thermal response of HVAC systems (e.g., [5]). EnergyPlus 
applies an explicit algorithm with non-coupled thermal and humidity bal¬ 
ances, and it assumes incompressible flow. The explicit algorithm uses a 
predictor-corrector approach to allow for large (typically 0.1 to 0.25 hr) time 
steps. The size of the variable time step is based on the limitation of the 
maximum temperature rise in any zone in the simulation, which ignores the 
physical time step limitation of air available for inhltration between adjacent 
zones. L-HVAC, on the other hand, couples the thermal transport, humidity, 
and density variations in flow in an implicit fashion, and its variable time step 
sizes are based on this inhltration limit. The coupling of humidity variation 
with thermal transport allows for direct incorporation of latent heat loads 
into the optimization process, therefore providing an advantage for L-HVAC 
over EnergyPlus for applications featuring high latent heat. The effects of 
coupling will be demonstrated later in this report. 

The L-HVAC algorithm comprises a set of layered nonlinear implicit so¬ 
lutions. All equations provided here are in SI units, but the code also allows 
for English units. Psychrometric calculations follow those provided in the 
ASHRAE Handbook of Fundamentals [16]. The code may be run in tran¬ 
sient or steady-state mode, and a transient run may be continued from a 
previous state. The time step used is variable and dehned as the minimum 
time for any lump in the simulation to transport all of its moist air to an ad¬ 
jacent downstream lump. Figure 3 graphically describes the following general 
algorithm: 
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• Solve the moist air flowrates and pressure drops through the duct sys¬ 
tem using a flow resistance calculation. These flow resistance calcula¬ 
tions include calculations of pressure drop, and temperature and hu¬ 
midity adjustments to the moist air flow due to the presence of the 

coil. 

• Solve the water mass balance in the system to obtain the humidity 
ratio distribution. 

• Solve the energy balance in the system to obtain the moist air temper¬ 
ature distribution. 

• Check for convergence; if not converged, then return to the flow resis¬ 
tance calculation with updated moist air density values. 

Each of these points are now described. 

2.1. Flow Resistance Calculations 

If the density of the moist air is approximated to be constant, then the 
system may be set up according to a flow resistance network calculation. 
This approach is an extension of that discussed in numerous undergraduate 
fluid mechanics texts (e.g. [17, 18]), and it allows for simplified analysis of a 
complex fluid network. In general, the pressure drop in a single duct branch 
(Ap) is given by 

Ap Q Rflow ^Pgain T ^Ploss (1) 

where Q is volumetric flow rate, Apgain is the pressure rise due to a fan, 
Apioss is a pressure drop due to the presence of a coil, and Rfiow is the flow 
resistance defined as 




1 (/ife + EA'.) 


( 2 ) 


where p is mass density, A is area, is hydraulic diameter, / is the Darcy 
friction factor, L is duct length, and Ki is a minor loss. In this equation, Ai 
and A 2 , refer to the area at the entrance and exit of the duct, respectively, 
and A is the average area of the duct. The Appendix provides details on 
the derivation of this equation. Note that the SI units for Rfiow are [kg/m’^j. 
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Figure 3: L-HVAC code algorithm. 
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which multiplied with volumetric flow rate [m^/s] squared yields a pressure 
drop in [Pa]. 

Coil calculations account for the Apioss term in Eq. 1. These calculations 
stem from reports by the Air-Conditioning and Refrigeration Institute (ARI) 
[19], and the equations account for the effects of partial or full condensation 
on the tubes, water loss due to condensation, pressure loss on both the air 
and refrigerant sides, and temperatnre rise in the refrigerant. Details of these 
calculations are too numerons to be inclnded in this article. The Apgain term 
stems from the application of a fan. 

The flow resistance calculation set includes pressure drop calculations 
along each branch following Eq. 1 and mass flow balance calcnlations at each 
lump: 




pinQin PoutQ out 


(3) 


where the subscripts in and out refer to the moist air streams into and ont 
of the Inmp, respectively. A general formulation for this setnp is written as 
the solntion for a system of nonlinear eqnations: 


duPi + di2P2 + ... + diMPM — RiiQi — R 12 Q 2 — ■ ■ ■ — RinQ'n ~ (4) 

Api 

OSS Apgain 

Riip\Q\ + D 12 P 2 Q 2 + • • • + Dinp'nQ'n ~ (3) 

^in ^in ^out Nout 

PiQiPjQ j 2 EE PiQiPjQ j 

i=l j=i+l i=l j=i+l 

where the solntion consists of a simultaneous solve of N instances of Eq. 4 for 
pressnre drop along a path and M instances of Eq. 5 for mass flow balance at a 
Inmp. The solution vector contains all M nnknown pressnres {pi,P 2 , ■ ■ ■ ,Pm) 
and N unknown flow rates sqnared {QijQl, ■ ■ ■ jQ'n)- The coefficients d 
and D in Eqs. 4 and 5 are either zero or one depending on the branch 
under consideration, and the coefficients R are either zero or the branch flow 
resistance Rfiow depending on the branch nnder consideration. The terms 
on the right-hand side of Eq. 5 acconnt for the nonlinearity of solving for 
flow rate sqnared when considering Eq. 3. Note that this adds a nonlinearity 
into the matrix solve, so a nonlinear solntion with a relaxation factor is 
implemented to allow for convergence. This convergence is achieved when 



the norm efiow falls below a user-defined value. The norm is dehned as 


€flo 


M 

E 

2=1 


iPi)k - iPi)k-i 


{Pi)k-1 


N 


2 = 1 


{Qi)k {Qi)k-1 


{Qi)k-l 


(6) 


where k is the iteration number. 


2.2. Water Mass Balance Calculations 

The flow resistance solution provides details regarding moist air flow rates 
in each duct in the system. Any changes in air humidity are accounted for 
in the water mass balance. In general, the increase in the amount of water 
niw in a lump is dehned by 

r\ ^in 

OTflyj , 

2 = 1 

where t is time, is the water mass how rate, rhinj is additional water 
injected into the lump, and and Nout are the number of moist air streams 
into and out of the lump, respectively. Equation 7 may be expressed in terms 
of humidity ratio ca, the mass of dry air in the lump rua, and the mass how 
rate of dry air ma as 


^out 

^ ^ ^ 121,2 "^inj (7) 

2=1 


d {maUj) 
dt 


^in ^out 

^ ^ kna,i^i ^ ^ rriiYij 

i=l i=l 


( 8 ) 


This equation, expressed in terms of psychrometric density p (kg of dry air 
per m^ of mixture), is 


d {pVoj) 

Ft 


^in ^out 

2=1 2=1 


(9) 


where V is the lump volume. Since the lump volume is constant, Eq. 9 can 
be restated as 

/^out \ / ^in \ 

( XjQij - ( j = (10) 
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where ^ = poo. For steady-state simulations or lumps with zero volume, this 
equation is written as 

/Nout \ /Ni„ \ 

^ j - ( j =^inj (11) 

The implicit form of the transient and steady-state equations is thus 

+ =(a+(w)(™-W 

V i=l / p+1 V i=l / p+1 ^ '' 

( 12 ) 

/Nout \ /Ni^ \ 

(Op+1 ( j ~ [Y1 ) = (^ini)p+i (13) 

\ i=l J p+1 \ i=^ / p+1 

for time step p with size At. A system of N equations is solved simultaneously 
to obtain all N unknown values of (Op+i- Note that values of (.^)p are 
updated at each iteration to account for density changes in the flow resistance 
calculation. 


2.3. Energy Balance Calculations 

The previous calculations provide values of moist air flow rates and water 
content in the system, which allows for determination of moist air properties 
in each lump based on psychrometrics and lump heating loads. A general 
energy balance in a lump follows 

0 7—1 ^in ^out 

f = (M) 

2 = 1 2=1 


where E is the energy in a lump, E is the energy flow into or out of a lump due 
to moist air flow, and q is the energy generated within a lump. The energy 
inside a lump may be expressed in terms of thermal and kinetic components: 

E = nia {C^a + OoCyw) T -\- rUa (1 -|- oj) 2 ^^ (^^) 

where T is temperature, C^a and C^yj are the air and water specihc heats 
at constant volume, respectively, and V is moist air speed. Dehning Cy = 
Cya + ^Cyy, and p = pCyT, Eq. 15 simplihes to 


E 




p + p{l + u)-V^ 


(16) 


10 



The energy flow into the lump E is 


E — rjQ + pQ (1 + Lo) 2 ^^ 


( 17 ) 


where Q = Q-pr and C„ = Cpa + ojCpu,, where Cpa and Cpw are the dry air 
and water speciflc heat at constant pressure, respectively. Combining Eqs. 
14, 16, and 17 for a fixed lump volume results in the expression 

Nir 


V— 

dt 


ri + p{l+uj)-V^ 


Z 


2=1 
r 


ViQi + PiQi (1 + aii) -V^ 


z 

2=1 


PiQi + PiQi (1 + -V^ 


+ q 


(18) 


The thermal terms are moved to the left side of the equation, and the kinetic 
energy terms are moved to the right side of the equation. Setting all 77 *, p*, 
and Ui for the exiting streams equal to p, p, and u, respectively, inside the 
lump yields 

^out ^in 


( P (1 + n;) 


d 


dt 


2=1 


2=1 


dt 


Ni, 


( PiQi (1 + ^i) l^^i 


Nout 


(19) 


- p (1 + ca) ( Qi-V^ ] + q 

\ - / j=i V / 

For steady-state simulations or lumps with zero volume, Eq. 19 is simplified 
to 

^out ^in ^in / -1 \ ^out 


P^Qi-Yl ^** 5 * = 2 

i=l i=l i=l ^ 


pi 


i=l 


(1 + ca) ( Qi-V^ ] + q 


( 20 ) 


The implicit formulation of Eqs. 19 and 20 at time step p is 

Ni„ 


Pp+1 


\ , ^out 


2=1 


p +1 


At 

V 


^ ^ PiQi 


= (P)» + 


2=1 


P+1 


P(l + + 2 ^^ 


P(l + + 2 ^^ 


p+i 

At 

V 


At 


'Ni, 


PiQi (1 + Y Yi 


. 2=1 


P+1 


^out 


pQ + ^)^Yi 


2=1 


p+1 


At. , , 

+ “^9p+i(21) 
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Vp+1 


^out 


i=l 


P+1 


Ni, 


^ ^ ViQi 


i=l 


'N„ 


'^PiQi (1 + (^i) 


P+1 


. 2 = 1 


P+1 


^out 


- p (1 + Cj) ^+ 9p+i (22) 


i=l 


P+1 


These iV equations are solved simultaneously for all N unknown values of 
? 7 p+i. As in the water mass balance calculations, the values of rjp are updated 
at each iteration to reflect changes in density and absolute humidity due to 
the flow resistance and water mass balance calculations, respectively. 


2.4- Convergence and Energy Consumption Calculations 

The Flow Resistance-Water Mass Balance-Energy Balance cycle continues 
until the solution converges below a user-dehned tolerance. L-HVAC dehnes 
the norm e at iteration k as 


e = 


M 

i=l 




Uu; 


'k-l 


+ 


M 

E 

2=1 


(+)a: (+)fc-l 




(23) 


where equals if > 0) 1.0 if = 0. 

L-HVAC is tailored toward HVAC systems that contain a chiller that 
distributes chilled water to various heat exchangers in the HVAC system. 
Based on this assumption, energy is input to three locations in the system: 


• Energy input to all pumps required to provide chilled water to the heat 
exchangers. 

• Energy input to all fans required to provide airflow. 

• Energy input to the chiller to cool the water. 


In L-HVAC, the pumps and fans are dehned based on one of 1) hxed pressure 
rise, 2) hxed howrate, or 3) a user-dehned fan curve. The pump/fan required 
power input is calculated as 


E. 


pumpi fan 


QAp 

^pump/ fan 


(24) 


where fjpumplfan 
lated as 


is the pump/fan efficiency. The chiller power input is calcu- 


E 


chiller — 


m^jCpnj (AT)^ 

COPchiller 


(25) 
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where COPchuier is the chiller Coefficient of Performance, and (AT)^ is the 
change in temperature of the chilled water as it passes through the chiller. 
Steady-state simulations provide values of energy usage per time and tran¬ 
sient simulations provide accumulated energy use. 

2.5. Controls 

The primary purpose of L-HVAC is to provide a means to predictively 
control HVAC conditions to allow for optimization per Figs. 1 and 2. There¬ 
fore, various types controls are available for implementation in the code. 

These controls include control dampers, variable fan/pump speeds, and vari¬ 
able chiller work input. The input mechanisms for these controls include 
thermostats and airflow regulators. Figure 4 shows the available pairing of 
input and output control mechanisms in L-HVAC. These controls follow a 
standard PID algorithm for an input voltage V based on deviation from 
setpoint At: 

V = Vo + +-^ J Aidi -|- q:£) ^ dt ^ 

where ap, «/, and ap are three user-dehned PID control coefficients. The 
use of PID algorithms is typical for HVAC modeling studies (e.g. [20]). Only 
ap is used in steady-state simulations. 

3. Validation Exercises 

Several tests were performed to ensure proper operation of the code, in¬ 
cluding a test comparing L-HVAC simulations of dry air to commercial code 
calculations and an incompressible theoretical model, a test showing behavior 
of control mechanisms, and a test demonstrating experimental-computational 
validation. Note that additional tests were used to show that the discretiza¬ 
tion of the ductwork and the time step had no signihcant effect on the sim¬ 
ulated results for the system described below. 

3.1. Theoretical Validation - Dry Air and Fixed Humidity Ratio 

L-HVAC was hrst validated by comparing its results to that by SINDA/FLUINT 
and a simplihed theoretical model for the system shown in Fig. 5. The sys¬ 
tem features 0.117 m^/s of ventilation air at 21.1°C mixed with 0.395 m^/s of 
return air flowing into a room of volume 28.3 m^. The room contains a load 
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of 1460 W. If density changes in the air are not taken into account due to 
pressure variation in the system, the predicted temperature rise in the room 
is governed by the differential equation 

dT 

pVCva-^ + xmCpaiT-To) = q (27) 

where Tq is the outside air temperature, and x is the volume fraction of 
external air into the room. The solution to this differential equation with 
T(0) = To is 


Tit) =To + 


q 

xmCpa 


1 — exp 


/ xmCpat \ 

V pC^aV ) 


(28) 


Figure 6 shows that all three approaches provide similar results for this simple 
system. The differences in the curves may be attributed to different applied 
approximations in the code algorithms or theoretical model. In general, the 
predicted response by L-HVAC agrees better with the theoretical model than 
the SINDA/FLUINT results. The maximum difference in predicted values 
between L-HVAC and the theoretical model and the values between L-HVAC 
and SINDA/FLUINT is 0.3°C and 1.3°C, respectively, which occurs at 3.5 
minutes. The average difference in values between L-HVAC and the theoret¬ 
ical model is 0.25°C, and the average difference in values between L-HVAC 
and SINDA/FLUINT is 0.48°C. Note that steady-state L-HVAC simulations 
obtain the same room temperature as transient simulations within O.Ol^C for 
t > 40 min. 

A modihcation of this test was used to ensure the uniformity of u when no 
water is injected or removed from the system. The same system of Fig. 5 was 
initialized with a globally uniform relative humidity of 50%, no mass injec¬ 
tion, and a 1460 W sensible heating load. The result shows that the system 
does not experience a change in uj greater than 0.002% of its initial value 
despite an increase in dry bulb temperature of 10.3"C, which is consistent 
with an overall mass balance of water in the system. 

Another interesting result is that u affects T. Analysis of Eq. 21 using 
Eq. 3 shows that even when kinetic energy terms are neglected, a lump ini¬ 
tialized with a globally uniform dry bulb temperature, no heating loads, and 
water mass injection will still experience a change in dry bulb temperature 
due to different airstream heat capacity values influenced by unique values 
of UJ in each airstream in the system. A modihed version of the dry air sim¬ 
ulation was produced to illustrate this fact. The sensible heat load in the 
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Figure 6: Theoretical validation via predictions of room temperature for a dry air system. 
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room was replaced with water injection of 0.756 kg/s. After 25 minutes, the 
room temperature dropped by 0.61°C, and the relative humidity increased 
to 59.4%. 


3.2. Code Validation to Experiment and Analytical Model 

A simple experiment and corresponding analytical model were used to 
validate the code. The lumped approximation is used for the analytical 
model. In the experiment, wet-bulb and dry-bulb measurements were taken 
during a classroom setting. The only air input to the classroom is through a 
wall-mounted fan to the ambient environment outside the building. Velocity 
measurements for the fan found an average input flow rate of 0.296 m^/s 
using a Tycos'^^vane anemometer, and an assumed ventilation rate of 0.5 air 
changes per hour per ASHRAE [21] provided an additional 0.0426 m^/s. The 
room volume was calculated from measurements to be 308 m^. The room 
exhaust was assumed to be the same as the input flow rate from the outside. 
An array of 16 ceiling fans were used to mix the air inside the room. Wet- 
and dry-bulb Temperature measurements were taken for 75 minutes using 
a Sierra-Misco™psychrometer, and psychrometric calculations were used to 
determine values of uj{t). 

The analytical model applies a known mass flow rate of air from outside 
the building into the classroom. The total and latent loads for each of the 
27 people in the room are 120 W and 44.0 W per person per ASHRAE [21]. 
The total amount of water injected into the room air is calculated from the 
latent heat to be 0.0291 kg/s using 


EfpQlat 

hfg 

where Np is the number of people in the room, niinj is the mass flowrate 
of water injected into the system, Qiat is latent heat, and hfg is the latent 
energy of water per mass. A simple differential equation for the average 
system humidity are expressed as 



pV 


doo 

dt 


m (wo - ca) ihinj 


(30) 


Eigure 7 shows that a reasonable agreement in humidity ratio variation is 
achieved if the external u is 0.0069, which corresponds to 75% relative humid¬ 
ity (RH) for 13°C air, which are reasonable external conditions for the time 
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of the experiment and considering the experimental error due to humidity 
gradients in the room and inhltration from other areas within the building. 
A dry bulb temperature validation is inadequate with this experiment as the 
effects from wall conduction, convection, and glazing effects provide too large 
of an experimental error. 

3.3. Application of Controls for Energy Savings 

Additional validation tests were performed to show the proper response 
by the simulated controls. In Fig. 5, a control damper is placed in the system 
at the inlet to the room. The room temperature setpoint is 26.7 °C for a 
cooling system containing a 3.517 kW chiller, a variable speed fan containing 
a single line fan curve connecting two points (700 Pa at 0 m^/s and 0 Pa 
at 0.7 m^/s), and a pump with 6.31 x 10“*^ m^/s output. The efficiency of 
the fan and pump is 0.8, and the COP for the chiller is 3.0. The external 
conditions are 21.2°C and 50% RH. Figure 8 shows that the system response 
to reach the setpoint depends on the output device. Each output device uses 
Eq. 26 (Up % 0 and «/ = ap = 0. 

Fig. 8 also shows that the steady-state temperature for a system without 
control is 23.6°C. A steady-state simulation of the system shows an energy 
use of 1.30 kW to attain this temperature. The energy savings per time for a 
room temperature setpoint of 26.7°C were calculated for system control via 
each mechanism using the formula 

E savings E^ic E(^ (31) 

where the subscripts nc and c refer to system power usage without controls 
and with controls, respectively. L-HVAC steady-state simulations can calcu¬ 
late energy savings for an inhnite variation of conditions, arrangements and 
control strategies. Results clearly show the dominant energy consumption of 
the chiller and the coupled nature of various control strategies. For example, 
varying the fan airflow will also reduce the air conditioning load since less 
reduced airflow over the cooling coil will reduce the amount of heat trans¬ 
ferred. A minimum quantity of airflow is required for proper circulation. The 
arrangement and cooling load variants will have a large effect on the energy 
savings for variable speed fans and some studies such as Teitel et al. [22] have 
reported savings as large as 25%. 
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Table 1: Effect of humidity on steady-state temperature prediction. 


Latent Heat Fraction 

Steady-State Temperature, °C 

0% 

31.5 

20% 

31.4 

40% 

31.3 

60% 

31.2 

80% 

31.2 

100% 

31.1 


4. Effect of Coupling Humidity and Temperature 

L-HVAC was designed to include the influence of humidity on tempera¬ 
ture, thereby allowing for optimization of systems featuring high latent loads. 
The effect of humidity on temperature is demonstrated by adjusting the frac¬ 
tion of latent to total heat generation in the room. Table 1 shows that a 
variation of 0.4"C in resultant steady-state room air temperature when the 
ratio of latent heat load to total heat load is increased from 0% to 100%. 
Although the uncertainty in approximating the room air as uniform super¬ 
sede this value, a hne discretization of a fluid system with high latent loads 
would beneht from incorporating the influence of humidity on the temper¬ 
ature calculation. This ability to incorporate the latent heat loads into the 
temperature calculation allows for the direct calculation of energy require¬ 
ments for both temperature and humidity control by the HVAC system, and 
therefore optimization of HVAC energy use by systems with latent loads is 
feasible using L-HVAC. 

5. Conclusion 

The L-HVAC code presented here provides the basis for optimization of 
HVAC designs to meet a particular load distribution. The incorporation 
of control mechanisms and energy calculations within the code allow for a 
direct means to observe the system energy usage for various control settings. 
The code organization allows for the user to specify the accuracy of the 
temperature, humidity, and energy calculations based on the discretization 
level of the model. Future work on this project will require hner-resolution 
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models than that shown in Fig. 5, but the general approach to modeling the 
system will not change. 
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7. Appendix: Derivation of Flow Resistance Relation 

The standard relation for incompressible, viscous internal flow between 
two axial locations in a duct (dehned here as locations 1 and 2) is [17, 18] 


P ^1/2 

-h a- -h ^ 1 — 

P9 ^9 , ^ 





Hl - Hgain + Hioss (32) 


where d, Hl, 9, and are the kinetic coefficient, duct head loss, gravitational 
acceleration, and elevation, respectively; and Hgain and R/oss are the addi¬ 
tional head gain and loss due to the presence of fans and coils, respectively. 
This equation may be rewritten as 


Ap = (pi - P2) = pg 


f/2 t/2 

0 : 2 -^ — + (^2 ~ Zi) + Hl 


^Pgain ^Ploss 


(33) 

where Apgain = P9Hgain and Apioss = P9Hioss- If turbulent flow is assumed, 
then di = 0(2 = 1- Furthermore, potential energy changes may also be 
assumed to be negligible due the low density of air. These assumptions 
simplify Eq. 33 to 


Ap = 


P 

2 



Apgain Apioss 


(34) 


where the head loss is expressed as 


Hl = 


‘29 



(35) 


Since Q = VA, Eq. 34 may be expressed as 


Ap = 




Aj A2 


D, 


Apgain Apio 


(36) 
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or 

Q Rflow ^Pgain ^Ploss 
where the flow resistance is deflned as 




1 1 

If ^ 13 



(37) 

(38) 
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Figure 7: Comparison of L-HVAC prediction of uj{t) to predictions by Eq. 30 and experi¬ 
mental measurements. 
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Time, min. 


Figure 8: Implementation of various control devices into a simulated HVAC system. Con¬ 
trols contain a setpoint of 26.7°C. The adjustable damper, variable chiller work input, and 
variable fanspeed contain values of ap of -1.0, 0.1, and 1.0, respectively. 
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